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Abstract 



We present submillimeter observations of various molecular ions toward two dense 
clouds, NGC 2264 IRSl and W 3 IRS5, in order to investigate their ionization frac- 
tion. Analysis of the line intensity ratios by the way of statistical equilibrium cal- 
culations allows determination of the physical parameters: n(H2) ~ (1 — 2) 10^ cm~"^ 
and Tkin ~ 50 — 100 K. Column densities and abundances are also derived. Together, 
the abundances of the observed ions provide a lower limit to the ionization fraction, 
which is (2 — 3) 10~^ in both clouds. In order to better constrain the electron abun- 
dance, a simple chemical model is built which calculates the steady state abundances 
of the major positive ions, using the observed abundances wherever available. With 
reasonable assumptions, good agreement within a factor of two with the observations 
can be achieved. The calculated electron fraction is Xe = (1.0 — 3.3) 10~® in the case 
of NGC 2264 and = (0.5 - 1.1) 10"^ for W 3 IRS5. In the first case, the high 
abundance of N2II"'" requires a rather high cosmic ray ionization rate > 10~^^ s~^, 
even if all nitrogen is assumed to be in gas phase N2. For W 3 IRS5, ionized metals 
such as Fe"'" and Mg"*" could provide 60% of the electrons. 
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1 Introduction 

The ionization fraction of dense clouds is a basic parameter of interstellar physics, 
since it determines the efficiency with which magnetic fields couple with the gas and 
ultimately controls the rate of collapse and star formation. Moreover, electrons and 
ions play a major role in interstellar chemistry, and their abundances provide infor- 
mation on the nature and importance of sources which ionize the gas. It is well known 
that in the diffuse gas, interstellar photons with A >912 A provide most of the ioniza- 
tion so that the electron fraction is of order 10""*, equal to the abundance of the major 
supplier C"*". This fraction is thought to drop by several orders of magnitude inside 
dense clouds, where the ionization mostly results from cosmic ray interactions. The 
major ions here are thought to be II3 , H"*", He"*", HCO"*", HsO"*" and possibly metal 
ions (e.g. Fe"*", Mg"*"). Unfortunately, except for HCO"*", all of these species are very 
difficult to observe directly for spectroscopic reasons. 

The first attempts to derive the electron abundance in dense clouds date back more 
than 15 years (Watson 1977; Guehn et al. 1977, 1982; Wootten et al. 1982). The basic 
idea was to use the observed DCO~'~/IICO~'~ ratio as a measure of the Il2D~'~/Il3 ratio. 
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which, together with a model for the deuterium fractionation and the H3 recombina- 
tion, gives a hmit on the electron abundance. Typical electron fractions derived in 
this way were 10~® — 5 10"'', with an order of magnitude variation from cloud to cloud 
(Langer 1985). However, this method was abandoned in 1984, when measurements 
by Adams & Smith (1984) found the dissociative recombination rate of H3 with e 
to be negligibly small, thereby invalidating the earlier determinations and raising the 
upper limits to 10~^. The H3 rate has subsequently been a heated topic of discussion 
in astrochemistry for more than a decade, preventing further studies of the ionization 
fraction. 

There are several reasons why a new attack on this problem is now both timely and 
warranted. First, the latest set of laboratory measurements on the H3 dissociative 
recombination rate seem to converge to the same (relative rapid) value under inter- 
stellar conditions (Amano 1990; Canosa et al. 1992; Larsson et al. 1993), although 
some lingering uncertainties remain (Smith & Spanel 1993). Second, many more ions 
can now be observed with the new sensitive, high frequency receivers: apart from 
HCO+ and its isotopes, these include HCS+, N2H+ and H3O+ (Phillips et al. 1992), 
whereas good upper limits on H2D"'" exist for a number of regions (van Dishoeck et al. 
1992). Third, the dissociative recombination rates of HCO"*", HCS"*", and N2H"'" have 
also recently been measured in the laboratory (Rowe et al. 1993). Fourth, the exci- 
tation and optical depth of the lines in dense clouds is now much better constrained 
by observations of several submillimeter lines, thereby minimizing the uncertainties in 
abundance determinations and physical conditions. 

We present here observations of various ions (HCO"*", H^'^CO"'", DCO"*", HCS"*" and 
N2H"'") in two different dense clouds, NGC 2264 and W 3 IRS5, in order to investi- 
gate their ionization fraction. The clouds are chosen to have a range in density and 
temperature (as determined from our earlier H2CO data) and to have complementary 
H30~'~ and/or H2D"'" emission observations available (Phillips et al. 1992; van Dishoeck 
et al. 1992), as well as infrared absorption line data on H3 (Black et al. 1990) and 
CO/H2 (Lacy et al. 1994). Previous studies in cold dark clouds such as TMC-1, with 
n(H2) = 10"* cm~"^, have found 5.8 10~^ 1^ Xe < 4.9 10~®, depending on the amount of 
PAHs and metals included in the modeling (Schilke et al. 1991). The sources studied 
in this work have higher temperatures and densities up to two orders of magnitude, 
and it will be interesting to investigate whether the electron fraction in them is lower, 
as predicted by theory. The observations and results are presented in Sections 2-3, 
and their analysis in Section 4. 

Together, the abundances of these ions provide a firm lower limit to the ionization 
fraction of the cloud. In order to constrain the electron abundance, more detailed 
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modeling is needed. The basic idea is that the molecular ions are formed by protona- 
tion of abundant neutral species in reactions with H3 : 

X + H+ ^ XH+ + H2 

where X=CO, H2O, N2, CS... Their destruction is governed by reactions with CO, 
H2O and and by dissociative recombination with an electron in the gas phase: 

XH+ + e ^ X + H 

Thus by measuring XH"*" and the abundances of the parents X for a sufficient number 
of species, it should be possible to constrain both the electron fraction and the H3 
abundance (or equivalently the cosmic ray ionization rate i^h) in dense clouds. This 
technique was pioneered by Wootten et al. (f979) based on the HCO~'~/CO ratio, but 
is extended here to include many more ion/neutral pairs. The actual modeling, which 
involves a more detailed network of reactions, is presented in Section 5. 

2 Observations 

The sources studied in this work are both associated with active regions of star for- 
mation. NGC 2264 IR is a well studied molecular cloud core [d ~ 800 pc) (Schwartz 
et al. f985; Kriigel et al f987). Its shape on an optical plate resembles a dark cone 
(the "Cone Nebula"). Presumably an intermediate mass star is forming here. The 
infrared source is bright (Allen 1972) and has been used as background light source 
for infrared absorption line studies of gas and solid phase species. In particular, an 
upper limit of 1.7 10~® exists on the II3 abundance in the gas (Black et al. 1990), and 
a lower limit of 1.7 10""* on the CO/II2 ratio (Lacy et al. 1994). 

W 3 IRS5 is the brightest infrared source [L = 1.7 10^ (Ladd et al. 1993)) in the 
W 3 Giant Molecular Cloud core [d ~ 2.3 kpc). A large amount of dense and warm 
material surrounds the massive star forming inside (Dickel 1980; Dickel et al. 1980; 
Hasegawa et al. 1994; Hayashi et al. 1989; Tieftrunk et al. 1995; Helmich et al. 1994). 
Because of the enormous visual extinction [Ay > 100 mag) the ultraviolet radiation 
from the young stars is not expected to influence the ionization balance deep inside 
the cloud. 

The observations of NGC 2264 were performed primarily at the James Clerk Maxwell 
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Telescope (JCMT^) in February 1994, with additional data taken in April, June and 
November 1994. The facihty receivers at 230 GHz (A2), 345 GHz (B3i) and 460 GHz 
(C2) were employed. These are double side-band receivers with upper and lower side- 
bands 3.00 GHz (A2, B3i) or 7.88 GHz (C2) apart. No systematic side-band checks 
were performed, since the spectra are simple and often there is only a single line in the 
spectrum. The image frequencies were compared against spectral-line catalogues and 
no obvious blending of lines was found. As the backend, the Digital Autocorrelation 
Spectrometer (DAS) was used in different spectral resolutions (mostly 0.094, 0.187 and 
0.374 MHz/channel). To improve the signal to noise, the spectra were smoothed to 
0.374 MHz/channel and in some cases to 0.750 MHz/channel. The former corresponds 
to ~0.5 km s~^ at 230 GHz. Integration times were generally 30 minutes on-|-off for 
230 and 345 GHz receivers but up to an hour for lines in the 460 GHz window, because 
of the higher system temperatures. This resulted in typical Icr noise in of 30 (A2)- 
50 (B3i) mK and 200-350 (C2) mK per resolution element (0.374 MHz). The pointing 
was checked regularly on OMC-1 and was found to be within 3". 

The calibration at the JCMT was performed with the chopper-wheel method. In 
February 1994, the beam efficiency was somewhat uncertain, because of inaccuracies 
in the shape of the telescope surface. However, careful comparison with earlier runs 
showed a reduction by only ~10%. This effect is corrected for by using an effective 
main beam-efficiency of 0.53, 0.45 and 0.35 for receiver A2, B3i and C2, respectively. 
During other runs, values of 0.72, 0.60 and 0.42, respectively, were used. In general 
we expect the absolute calibration to be accurate to about 30%; relative uncertainties 
can be significantly smaller for all runs. The beam size of the JCMT at 230, 345 and 
460 GHz is ~20", 15" and 12", respectively. For most spectra, position-switching by 
600" was found to be appropriate, except for the ^^CO lines for which a larger switch 
was used. Additional data on other species toward NGC 2264 IRSl will be presented 
by Schreyer et al. (1996). 

A limited number of spectra of NGC 2264 were taken at the Caltech Submillimeter 
Observatory (CSO^) in January 1994 with the 230 GHz (H2CO at 218 GHz) and 
345 GHz (C^^O 3-2) receivers. The low (500 MHz) and high (50 MHz) resolution 
acousto-optical spectrometers were used as backends. Pointing is reproducible within 
5" on the CSO. Main beam efficiencies of 0.72 and 0.60 were adopted, respectively. 
The CSO data at 218 GHz refer to a beam of ~32"; those at 345 GHz to a 20" beam. 

^The James Clerk Maxwell Telescope is operated by the Observatories on behalf of the Particle 
Physics and Astronomy Research Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 

^The Caltech Submillimeter Observatory is operated by the Californian Institute of Technology 
under funding from the U.S. National Science Foundation (AST 93-13929). 
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Earlier observations by van Dishoeck et al. (1992) with the CSO have resulted in 
upper limits on the H2D"'" abundance. These same spectra also reveal surprisingly 
strong N2H"'" 4-3 emission in this source in a 18" beam. 



NGC 2264 
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The observations presented here for W 3 IRS5 form part of the large JCMT 345 GHz 
survey of three star-forming cores in the W 3 cloud (IRS4, IRS5 and W 3(H20)) 
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Figure 1: Observed spectra toward: a and b NGC 2264 IRSl and c W 3 IRS5. Note 
that the DCO"*" 3-2 hne toward NGC 2264 consists of two components. Only the 
broad component was used in the analysis. The N2H"'" 5-4 line represents the first 
detection of this line in the interstellar medium. The N2D"'" 3-2 line may be present at 
the 2(7 level. Since W 3 IRS5 is a well-known outfiow source, only the central Gaussian 
was used in the analysis. The line next to the HCS"*" 8-7 feature at Vlsr = —20km s~^ 
is the ■^'*S02 19i^i9 — 18o,i8 line from the image side-band. 
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between January 1992 and November 1994. A complete log and details of the obser- 
vations will be given in a separate publication concerning the spectral survey (Helmich 
& van Dishoeck 1995; see also Helmich et al. 1994 for first results). For all runs, the 
facility receivers (Al, A2, B3i and C2) were used. Prior to July 1993, the acousto- 
optical spectrometer AOSC (500 MHz) was employed as the backend; subsequently 
the DAS was used in different spectral resolutions, mostly in its 500 MHz bandwidth 
mode. Binning of two channels resulted in a typical resolution of 0.750 MHz/channel. 
As for NGC 2264, no systematic side band identification was performed, because the 
line crowding in IRS5, compared with the two other W 3 sources in the spectral-line 
survey, is such that there is little ambiguity and possible blends are easily identified. 
The pointing was checked regularly on the nearby continuum source W 3(0H) and is 
within 3". The observations of W 3 IRS5 were performed in the 180" beam-switch 
mode, except for the ^^CO lines where a larger position switch was used. 

Additional observations of the HsO"*" ion toward W 3 IRS5 were obtained by Phillips 
et al. (1992) using the CSO, whereas an upper limit on H2D~'~ was obtained by van 
Dishoeck et al. (1992). The latter study also reported detection of the N2H"'" 4-3 line. 

3 Results 

Tables 1 and 2 summarize the resulting Gaussian fit parameters for the various lines 
for the two sources. Examples of spectra are given in Fig. 1 a, b and c. Although the 
profiles of the lower lines of DCO"*" and other species toward NGC 2264 suggest the 
presence of two components (Guelin et al. 1982; Schreyer et al. 1996), it was found 
that the lines observed in this work are well represented by single Gaussians with little 
spread in the line widths, the optically thick ^^CO 4-3 and HCO"*" 4-3 lines being 
exceptions. Most other lines are optically thin, or at most moderately optically thick. 
For example, the C'^'^S 5-4 line is weak, and the intensity ratio CS/C"^'*S=16-21 is only 
slightly smaller than the cosmic p^S]/p'*S] abundance ratio of 22.59. As mentioned 
above, N2H"'" is very strong in this source; the 5-4 line at 465 GHz presented in Fig. lb 
is the first detection of this line in any astrophysical source. N2D"'" 3-2 appears present 
as well, but only at the 2cr level. 

The lines toward W 3 IRS5 are also well represented by single Gaussians but there is a 
much larger spread in the line widths. The molecules HCS"*", C"^"^S and C'^'^S have much 
narrower lines than e.g., H^'^CO"'", HCO"*" and CO-isotopomers. It is easy to explain 
the larger line widths for CO and HCO"*" by optical depth effects, but this becomes 
more difficult for the rarer species. Part of the explanation may lie in inaccurate 
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Table 1: Parameters ol Gaussian fits lor NGC 2264 



Line 


Tmb 


Ay 


/ TuBdV 




(K) 


km s~^ 


K km s~^ 


CO 4 — 3 


54.3 


10.1 


586 


CO 2 — 1 


6.74 


3.0 


21.1 


C^'O 3 — 2 


4.18 


2.3 


10.3 


CS 5 — 4 


8.48 




3.3 


30 


CS 7 - 6 


5.42 


3.1 


18 


CS 10 — 9 


< 2.44 






C-^^S 5-4 


0.53 


2.5 


1.42 


HDO loi - Ooo 


< 1.94 






HCO+ 4 — 3 


13.6 


4.1 


59.3 


H^^CO+ 3 — 2 


1.34 


2.8 


4.0 


H^^CO^ 4 — 3 


0.96 


2.1 


2.1 


r)CO+ 3-2 


0.42 


3.2 


1.4 


HCS+ 6-5 


0.26 


2.2 


0.64 


N2H+ 5-4 


4.1 


2.5 


10.9 


N2D+ 3-2 


0.11 


0.8 


0.08 


[C l]3Pi -3 Po 


20.6 


4.5 


98.9 



Coordinates (B1950.0): a = 06 38 25.0, (5 = 09 32 29. 
All data were taken with the JCMT except lor C^^O 3-2. 
The upper hmits are 2cr in a 0.625 MHz channeL 

determinations ol the C"^"^S and HCS"*" fine widths, because these fines are so weak. 
However, source structure can play a role as well. A complete investigation is beyond 
tfie scope ol tfiis paper; tfierelore, lor eacfi species its observed line widtfi is adopted 
in tfie analysis. Tfie detection ol C"^"^S implies tfiat even tfie C'^'^S lines are sligfitly 
optically tfiick. From tfie ^"^CS line, fiowever, it can be seen tfiat tfie optical deptfi 
must be small. 
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Table 2: Parameters ol Gaussian fits lor W 3 IRS5 



Line 




Tmb 
(K) 


Ay 

km s~^ 


J iVlJ_> 

K km s~^ 


CO 3-2 




77.7 


13.2 


1099 


CO 2-1 




103 


10.7 


1177 


C^^O 2 - 


- 1 


11.2 


4.59 


54.7 


C^^O 3 - 


- 2 


3.87 


4.09 


16.9 


CS 7 - 6 




6.53 


3.66 


25.4 


C^^S 5 - 


4 


0.19 


2.84 


0.58 


C^^S 7 - 


6 


0.59 


1.35 


0.85 


C-^^S 5 - 


4 


1.14 


2.91 


3.54 


C-^^S 7 - 


6 


0.77 


1.83 


1.5 


^^CS 5 - 


4 


0.35 


2.77 


1.04 


HCO+ 4 


_ 'X 


24 8 


6 1 ^ 


1 62 


H13C0+ 


3-2 


2.06 


4.21 


9.28 


H13C0+ 


4-3 


2.49 


3.47 


9.18 


HCi^O+ 


4-3 


0.15 


3.3 


0.56 


DCO+ 5 


-4 


0.22 


2.67 


0.64 


HCS+ 6 


- 5 


0.23 


0.88 


0.21 


HCS+ 8- 


- 7 


0.28 


1.58 


0.46 


N2H+ 5 - 


- 4 


<0.38 






N2D+ 3 - 


- 2 


< 0.1 






[C i]3Pi 


3 p 


16.5 


7.05 


124 



Coordinates (B1950.0): a = 02 21 53.1, S = 61 52 20. 
All data were taken with the JCMT. 
The upper hmits are 2cr in a 0.625 MHz channeL 
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4 Analysis 



4.1 Physical parameters 

In order to derive the column densities from the observed line intensities, first the 
physical parameters in the sources have to be determined. Statistical equilibrium 
calculations were performed along the lines described by Jansen et al. (1994) and 
Jansen (1995). These models take the various coUisional and radiative excitation and 
de-excitation processes into account, and treat the radiative transfer with an escape 
probability method. Because the observed lines are mostly optically thin, the precise 
details of the radiative transfer method do not affect the results. The ratios of some 
H2CO lines are particularly sensitive to temperature variations (see e.g. Jansen et al. 
1993; Helmich et al. 1994; Mangum & Wootten 1993) and are thus very well suited 
for kinetic temperature determinations. Lines from linear rotors such as CS, C'^'^S, or 
ions like H^'^CO"'", on the other hand, are more sensitive to density. 

Line intensity ratios were calculated for a range of densities and temperatures for 
H2CO, CS, H^'^CO"'" and N2H"'", and are plotted in Fig. 2a, together with the observed 
ratios for NGC 2264 (see also Schreyer et al. 1996). The contour plot of the CS 7-6/5- 
4 line ratio is produced for the slightly optically thick case appropriate for this source. 
All other ratios refer to optically thin lines. From the H2CO 3o3 — 2o2 and 322 — 221 lines 
at 218 GHz measured in the same frequency setting, a kinetic temperature of 55^5° K is 
found. Furthermore it is seen from Fig. 2a that the measurements for the N2H"'" and CS 
ratios indicate a similar density of (2 — 3) 10^ cm~"^. Thus, these species probably probe 
the same region. The density indicated by H^'^CO"'" appears to be somewhat lower 
than that found from the two other species. This difference is difficult to account for 
in terms of optical thickness, but could be due, e.g., to uncertainties in the coUisional 
rate coefficients or source structure. The measurements are, however, within the errors 
in agreement with each other and give a number density of n(H2) = (2 it 1.5) 10^ cm~"^. 
Maps of the CS 7-6 and 5-4 and ^^CO 3-2 lines by Schreyer et al. (1996) show that 
the temperature cannot be much lower than that found from H2CO. Note that the 
above single temperature and density analysis assumes that all species and lines are 
similarly distributed in the beam. No corrections are made for the fact that the higher 
frequency JCMT data refer to a smaller beam than those at lower frequencies. 

For W 3 IRS5 the temperature Tkin ~ 100^20 K and density n{E2)^ 10^ cm ^ derived 
in Helmich et al. (1994) from various formaldehyde lines was adopted. The H^'^CO"'" 
4-3/3-2 and C'^'^S 7-6/5-4 line ratios, presented in Fig. 2b, show that these molecules 
probe about the same density, suggesting that they are present in the same gas. It is 
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Figure 2: Contour plot showing ratios ol lines. Shaded regions refer to the observed 
range toward the source: a NGC 2264, b W 3 IRS5. All ratios refer to optically thin 
lines, except the CS 7-6/5-4 ratio toward NGC 2264, where iV(CS) = 1.3 10^^ cm'^ 



and Ay = 3.2 km 



was 



adopted. The N2H"'" ratio for W 3 IRS5 is uncertain and 



is therefore plotted as an upper limit. 
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assumed that the emission of all molecules described here fills our beams. 

The inferred densities toward our two sources are indeed two orders of magnitude 
larger than those found in dark clouds such as TMC-1, which have been the subject 
of previous studies of the electron fraction. 

4.2 Column densities 

The total beam-averaged column densities for each species are calculated by fitting 
the absolute intensities, using the temperature and density derived from our statistical 
equilibrium calculations. The results are given in Tables 3 and 4 for the two sources. 
Since the HCO"*" 4-3 line is optically thick toward NGC 2264, the column density 
of HCO"*" is determined from H^'^CO"'", using the [^^C/^"^C] ratio of 94 obtained by 
Mitchell & Maillard (1993) for this source. The ^^CO column density was calculated 
from the C^''0 3-2 line assuming that the CO/C^''0 ratio is the same as the cosmic 
[^^0]/[^''0] ratio of 2600. The lower limit on the CO/H2 ratio found by Lacy et al. 
(1994) of 1.7 lO""* then gives an upper limit on the total beam-averaged column- 
density A^(H2)~ 8.3 10^^ cm~^. Although the CO/H2 ratio derived from the pencil- 
beam infrared absorption line observations does not necessarily apply to the whole 
volume sampled by the emission data, the CO column density derived from the two 
methods agree within a factor of two. 

The total (beam-averaged) column density toward W 3 IRS5 is derived from the ''O 
3-2 line. In contrast with Helmich et al. (1994), a CO/H2 abundance ratio of 2.7 10""* 
was used instead of 8 10~^. The value in Helmich et al. (1994) was chosen because the 
study of Tielens et al. (1991) of the solid CO seemed to indicate a lower value. The 
subsequent study of Lacy et al. (1994) however suggests a much larger ratio for warm 
sources like NGC 2024 IRS2. Since the gas and grain temperatures toward W 3 IRS5 
are even higher than those toward NGC 2024, a higher value is probably more appro- 
priate. This value decreases the total beam-averaged column density from 5 10^"^ cm~^ 
to 1.3 10^"^ cm~^, and at the same time increases the abundances by a similar factor. 
It is stressed again, however, that the derived column-densities and abundances are 
uncertain for this source due to the unknown structure and distribution of the various 
species within the beam. The relative values should be more reliable, however. 

For both clouds the observation of the [C l] 492 GHz line provides good upper limits 
on the beam averaged column density and abundance of atomic carbon since a large 
fraction of this emission probably arises from the outer skin of the cloud. 

Tables 3 and 4 include column densities and abundances (w.r.t. H2) for other relevant 
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Table 3: Column densities and abundances in NGC 2264 



c 

opecies 


A T 

iV 


X 




(cm j 




CS 


1.3 10^4 


1.6 10-9 


c^^s 


4.0 10^2 


4.8 10"^^ 


HDO 


< 7.1 10^3 


< 8.6 10"^° 


H13C0+ 


1.6 10^2 


1.9 10"^^ 


HCO+ 


1.5 


1.8 10-9 


DCO+ 


7.4 10^^ 


8.9 10-^2 


HCS+ 


1.8 10^2 


2.2 10"^^ 


N2H+ 


2.5 10^3 


3.0 10"^° 


N2D+ 


4.6 10^° 


5.5 10-^3 


H2D+ 


< 9.0 10^2 


< 1.0 10"^° 


H3+ 


< 1.4 10^^ 


< 1.6 10"^ 


H2O 


< 4.2 10^6 


< 5.1 10"^ 


c 


< 1.8 10^^ 


< 2.2 10"^ 



species obtained from the literature for the two sources. In the case of NGC 2264 
we adopted the H2D~'~ column density derived by van Dishoeck et aL (1992), with 
shghtly different total H2 density and temperature. The H2O column density is taken 
from Wannier et al. (1991) who derived H2O/CO < 0.003 in NGC 2264 from searches 
for H2^^0 using the Kuiper Airborne Observatory in a 2.7' beam. For HsO"*" toward 
W 3 IRS5, the data given in Phillips et al. (1992) were reanalysed with the above 
density and temperature. From observations of the H2O 183 GHz line presented by 
Cernicharo et al. (1990), Phillips et al. (1992) determined a;(H20) ~ 10"^ - 10"^ The 
O2 limit of < 10^9 cm~^ for that source is based on a deep search for the ^^0^*0 234 
GHz line (Keene, van Dishoeck & Phillips, unpublished result). 

5 Models 

In order to estimate the electron fraction, we built a simple chemical model which 
calculates the steady state abundances of the major positive ions, namely: H"*", H^, 
H+, He+, C+, H3O+, HCO+, DCO+, N2H+, H2D+, HCS+, N2D+ and H2DO+. The 
electron density is set equal to the sum of the densities of the positive ions, including 
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Table 4: Column densities and abundances in W 3 IRS5 



c 

opecies 


AT 
iV 


X 




(cm j 




cs 


2.0 


1.5 10-9 




1.4 10^^ 


1.1 10-^ 




3.3 10^2 


2.5 10-11 


HCO+ 


2.0 


1.5 10-9 


DCO+ 


4.8 


3.7 10-12 


HCS+ 


6.3 


4.8 10-12 


N2H+ 


< 8.3 


< 6.4 10-12 


N2D+ 


< 3.5 


< 2.3 10-12 


H3O+ 


5.6 10^3 


4.3 10-1° 


H2D+ 


< 5.8 10^2 


< 4.4 10-11 


O2 


< 1.0 10^^ 


< 7.7 10-^ 


c 


< 2.1 10^^ 


< 1.6 10-^ 



metals or PAH ions. Only gas phase reactions are included. The temperature, density 
of molecular hydrogen and abundances of CO and CS are those derived from the 
observations for each cloud. 

The abundances of 0, O2, H2O, N, N2, D and HD are free parameters, taking into 
account the constraints imposed by other observations or models, and allowing for 
a certain amount of depletion on grains. The main motivation of this approach is 
that detailed chemical models generally give O2 and H2O abundances higher than 
the upper limits derived from the observations [x ~ 10-^ — 10-^) (see Tables 3 and 
4). Moreover these models usually do not allow for grain surface reactions and do 
not consider deuterated varieties. In our case the abundance of O2 does not play 
an important role since this species contributes solely to the destruction of He"*", H"*" 
and C"*", for which we do not have observed abundances, and for which the calculated 
abundances are low in any case. Other free parameters are the cosmic ray ionization 
rate i^h and the abundance of ionized metal Xy[+ . 

Concerning the reaction rates, special attention has been paid to the dissociative re- 
combination rate coefficients of HCS"*", HCO"*", N2H"'" and H3 , which have recently been 
studied in the lab. For HCO"*" and H3 , the rate coefficients measured by Rowe et al. 
(1993) and Amano (1990) are very close: from their work, we derive the temperature- 
dependent values: 
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A;hco+ = 2.9 10-^(r/300)-°-^^ cin^s-\ 

k^^ = 1.8 10-'(r/300)-°-^^ cm^s-^ 

In fact, Smith & Spanel (1993) obtain a lower value = 3 10~®cm"^s ^ at 300 K, 
but use of this value does not affect our results since H3 is mainly destroyed by CO, 
due to the high density and the high abundance of CO. 

For N2H"'", the rate measured by Amano (1990) can be fitted by: 

A;n^h+ = 7.3 10-^(r/300)-°-'^ cm^s-^ 

For HCS+, Rowe et al. (1993) found a rate coeflicient of (6 - 8) 10"^ cm^s'^ at 300 
K, which is about one order of magnitude larger than found in previous studies (Millar 
et al. 1985; Talbi et al. 1989). The original motivation for slow HCS"*" dissociative 
recombination came from the observed high HCS~'~/CS abundance ratios. Rowe et al. 
(1993) claim that they can still reproduce this high value with a rapid rate by assuming 
that the branching ratio toward the H + CS channel represents only one third of the 
total dissociative recombination process. Since, in addition to this assumption, the 
infiuence of the vibrational excitation of the HCS"*" ions in their experiment is unclear, 
the original determination of Millar et al. is kept: 

A;hcs+ = 5 10-^(r/300)-°-' cm^s"^ 

For the deuterated species (DCO"*", H2D"'", N2D"'" and H2D0~'~) the same dissociative 
recombination rate coefficients as for the hydrogen-based ones have been adopted; the 
main formation and destruction rate coefficients are those given by the UMIST group 
(Millar et al. 1991). In particular, we followed the study of Millar et al. (1989) for the 
temperature-dependence of the important exchange reaction, II3 -|- HD — )• II2D"'" -|- II2: 

ki = 1.8 10-^(1. - r/340) cm^s"^ 

For the reverse reaction, Pagani et al. (1992) have shown that the rate depends on the 
ortho-Il2D"'"/para-Il2D"'" ratio. Since it is beyond the scope of this study to calculate 
this ratio in our model, we simply adopt: k^ = k^e'^^'^^'^ . 
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The best fitting models for the two clouds are described in the next section. The 
model parameters are summarized in Tables 5 and 7. 



The principal characteristics of the NGC 2264 observations are the high abundance 
of N2H"'", and the high HCS~'~/CS ratio equal to 1.4 10~^. As noted by van Dishoeck 
et al. (1992), the high N2H"'" abundance implies that virtually all nitrogen must be in 
the gas phase in the form of N2. Also, the CO/H2 ratio has to be close to the lower 
limit derived by Lacy et al. (1994), i.e., CO/H2 = 1.7 10"'*, because the dominant 
destruction route for N2H"'" is by reaction with CO, and the main formation route 
is by: N2 + H3 — )• N2H"'" + H2, and H3 is also mainly destroyed by CO. H3 , in 
turn, is formed through cosmic ray ionization of H2, and its abundance is directly 
proportional to i^h (Lepp et al. 1987; Black et al. 1990; Phillips et al. 1992). 

As can be seen from Tables 5 and 6, the best fit requires that all nitrogen is in gas 
phase N2 for a rather high value of i^h ~ 2.5 10~*^ s~* and low abundances of atomic 
oxygen and carbon (Model 1). In fact, if x{0) ^ 8.5 lO""* and x{C) ^ 1.4 lO""*, the 
abundance of N2H"'" can be written as: 



Thus, the N2H"'" abundance can also be reproduced by an even higher value of i^h, but 
with the fraction of nitrogen in N2 decreased. The dependence of the ion abundances 
on (u is illustrated in Fig. 3a. The inferred values for i^h are larger than those normally 
assumed for dense clouds, i^h ~ (1~5) 10"*'' s~*, but could be due to e.g., the presence 
of X-rays from the young stellar object, which have a similar effect on the ionization 
(Casanova et al. 1995, Farquhar et al. 1994). If the abundances of and C are close 
to the upper limits derived from observations and other estimates (Little et al. 1994, 
1995; Pagani & Fncrenaz 1994) i^h has to be very high: an extreme case is given by 
Model 2, in which we have assumed x[0) = x[CO) and x[C) = 0.1 x[CO). 

The main problem comes from reproducing the high HCS"*" abundance. It can be 
shown that its abundance can be expressed as: 



5.1 NGC 2264 IRSl 
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Figure 3: Ratios of calculated ion abundances to observed ones as functions of the 
three dominant parameters: i^h, x(H20) and x(M"'"): a in the case of NGC 2264 (Model 
1), b in the case of W 3 IRS5 (Model 2). 
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Table 5: Model inputs for NGC 2264 



r Tri m pf pr 


Model 1 


Model 2 


Ch (s-^) 


2.5 10-^6 


3. 10-45 


x{M+) 


< 2 10"^ 


< 10-^ 




1.1 10-4 


3.1 10-5 


X[tl2U ) 


iU 


iU 


x{0) 


< 5. 10-6 


1.7 10-4 


x{C) 


< 5. 10-6 


1.7 10-5 


x{B) 


2.8 10-^ 


8. 10-^ 


x(HD) 


2.8 10-5 


10-5 



so that a;(HCS~'~) increases when the electron fraction decreases: it is one of the rea- 
sons why the metal ion abundance has to be low in this case. Moreover, keeping in 
mind that the abundance of H3 is now fixed - because i^h is fixed from the N2H"'" 
observations - , a;(HCS~'~) depends essentially on the abundances of HCO"*" and HsO"*", 
where a;(HCO~'~) has been constrained from observations. As can be seen in Fig. 3a, 
the dependence of the HCS"*" abundance on i^h is not large, because in the model HsO"*" 
is the dominant ion leading to the formation of HCS"*". 

The formation of HCO"*" invokes CO, H3 and N2H"'" . Its destruction proceeds mainly 
by dissociative recombination and by reaction with H2O and atomic deuterium, which 
results in DCO"*". HsO"*" is formed by reactions of H2O and with molecular ions 
(essentially HCO"*" and H3 ) and is destroyed by dissociative recombination. 

At this point, the important remaining parameters are thus the abundances of H2O, 
D and HD which govern the abundances of HCO"*", HsO"*" and DCO"*", the latter 
via HCO"*" and H2D"'". Fig. 3a shows the model HCO"*" and DCO"*" abundances as 
functions of the H2O abundance. In order to reproduce the observed values a rather 
high abundance of II2O of at least 3 10-5 is needed. Also, it is inferred that only 1 to 
5% of deuterium is in atomic form. The high value of a;(Il20) seems in contradiction 
with the work of Wannier et al. (1991). However, their results refer to a much larger 
beam of ~ 2.7' than the region sampled by our observations. The presence of water 
masers near IRSl suggests that locally the II2O abundance may be significantly higher 
(Henning et al. 1992). 

Finally, as can be seen in Table 6, good agreement between the observed and calculated 
abundances is achieved for N2H+, HCO+, DCO+, N2D+ and HCS+. 
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Table 6: Model fits for NGC 2264 



x (model) /xfobs) 


Model 1 


Model 2 


HCO+ 


1 


1 7 


DCO+ 


0.8 


1.1 


N2H+ 


1.0 


1.0 


N2D+ 


1.4 


2.4 


HCS+ 


0.7 


0.8 


predictions 


X Q 


9.6 10-9 


3.3 10-^ 


x{Et) 


3.9 10-1° 


2.9 10-9 


x(H2D+) 


6.7 10-13 


1.9 10-12 


a:(H30+) 


6.8 10-9 


2.5 10-^ 


x(H2D0+) 


2.0 10-11 


3.1 10-11 


a:(NH3) 


5.1 10-9 


3.1 10-9 


x(NH2D) 


1.9 10-12 


5.0 10-13 



With these assumptions, the electron fraction is Xe = (1.0 — 3.3) 10-®, the dominant 
ion being H30~'~. Searches for this ion in this object would therefore be very valuable. 
A strict lower limit is provided by the sum of all observed ions: Xe > 2.2 10-9. The 
maximum value of a;(M"'"), for which the ratio of the calculated to observed abundances 
lies between 0.5 and 2, is 4 10-9 in the case of Model 1 (Fig. 3a) and 2 10~^ in the case 
of Model 2. The predicted abundance of H3 is below the upper limit derived by Black 
et al. (1990) and is in good agreement with the values estimated by van Dishoeck et 
al. (1992) from H2D"'", in their high density case, especially when it is assumed that 
all nitrogen is in gas phase N2. Model 1 predicts H2D"'"/H3 and H2DO~'~/H30~'~ ratios 
comparable to the observed ratios DCO+/HCO+ ~ 5 10-3 and N2D+/N2H+ ~ 2 10-3 
(Table 9). In the case of Model 2, the calculated N2D"'"/N2H"'" is higher due to the 
higher abundance of D, and the H2D"'"/H3 is lower due to the lower abundance of HD. 
The model also predicts the abundance of NH3, which is close to the value of 3 10-9 
determined observationally by Schreyer et al. (1995), and that of NH2D. It should 
be noted, however, that the current model does not include grain surface reactions, 
which may be significant for the production of NII3 and its deuteration. 
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5.2 W 3 IRS5 



The excitation analysis ol §4 has shown that the density in W 3 IRS5 is slightly 
lower than that in NGC 2264, but the temperature higher. Moreover the observed 
HCS~'~/CS ratio of 3 10""^ and the N2H"'" abundance are lower, suggesting a higher 
electron fraction and/or a lower cosmic ray ionization rate. The adopted CO/H2 
ratio is also higher than toward NGC 2264. The model analysis for W 3 IRS5 cannot 
proceed in the same way as for NGC 2264, because no good determination of the N2H"'" 
abundance is available. Although N2H"'" 4-3 is most likely detected by van Dishoeck 
et al. (1992), the line falls at the edge of the spectrum and does not lead to a reliable 
column density. Moreover, the 5-4 line was not detected. On the other hand, HsO"*" 
data are available for this source, in contrast with NGC 2264, which together with 
HCS"*" proved to be very useful in the analysis. 

The abundance of HCS"*" can be written as: 



Using the observational constraints on HCS"*", HsO"*" and HCO"*", we can derive a 
relation between ^(Hg ), x{0) and Xg from this expression. 

The relation between Xg, x{}i^) and x{M~^) is given by: 



Together, these two formulae lead to an expression of ^(Hg ) as a function of x{0 ) 
and x{M+). 

Further relations can be derived from 




xe ~ x(H30+) + x(HCO+) + x(H+) + x(M+), 



or 



~ 2 10"^ + x(H+) + x{M+). 





[x(H20)(x(HC0+) + 1.7x(H+)) 
+0.3x(O)x(H+)] 
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which leads to a;(H20) as a function of x[0 ) and x(M'^). In the same way an 
expression for x(D) as a function of x[0 ) and x(M'^) is derived from: 



x(N2H+) + 1.9x(H+) 

I ^(H20) I xj-D) 
9.7 10-9 T 2.5 10-6 T 5 3 ;^q-6 

Assuming that a;(HD) = [D]/[H2]— x(D), with [D]/[H2] = 2.8 10~^, another expression 
of x(D) can be obtained from: 



x(HCO+) 
1.5 10-9 



x(DCO- 



3.7 10-12 



227.3 



^(H20) 



:d)(x(hco+) + 



7.7x(H+)) + 10x(HD)x(H+)] 



The two expressions of x(D) imply that: 



x(0) 



3.9 10- 



+ 



x{M+) 
8.2 10-9 



1. 



The maximum value of x(M'^) is then around 8. 10-^, corresponding to a low atomic 
oxygen abundance (Model 2). The values of these two parameters being fixed, we 
can obtain the abundances of H2O, D and HD needed to reproduce the observations, 
as presented in Table 7: the abundance of H2O adopted in the modeling is in good 
agreement with the value derived by Phillips et al. (1992) from data of Cernicharo et 
al. (1990), and only 0.04 to 0.1% of deuterium is in atomic form. We also derive the 
abundance of H3 from x{0 ) and x{M~^), and then i^h since: 



2.1 IO^Ch 



1 + 



5.7 10-* 



+ 



2.3 10- 



+ 



9.6 10-^ 



In the case of Model 1 the values of x[0) and x(M'^) have been chosen in order to get 
a value of i^h close to the lower limit generally assumed in dark clouds. 

As for NGC 2264, the sensitivity of the calculated ion abundances to the main pa- 
rameters, a;(M"'"), a;(Il20) and i^h are shown in Fig. 3b, in the case of Model 2. The 
fraction of N2 is chosen to give an abundance of N2II"'" close to the observed upper 
limit: between 20 and 40% of the nitrogen is then in gas phase N2, similar to the 
values suggested by van Dishoeck et al. (1992). 

With the values of the parameters (Table 7), good agreement between the observed and 
calculated abundances is achieved for HCO"*", DCO"*", N2ll~'~, HCS"*" and HsO"*" (Table 
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Table 7: Model inputs for W 3 IRS5 



Parameter 


Model 1 


Model 2 


Ch (s-^) 


10-1^ 


2.5 10-^^ 


x{M+) 


3 10-9 


8 10-9 




4.4 10-5 


2 10-5 


x(H20) 


10-6 


3 10-6 


x{0) 


2.5 10-5 


< 2. 10-6 


x{C) 


< 2. 10-6 


< 2. 10-6 


x{B) 


1.2 10-^ 


3 10-^ 


x(HD) 


2.8 10-5 


2.8 10-5 



Table 8: Model fits for W 3 IRS5 



a;(model) /a;(obs) 


Model 


1 


Model 


2 


HCO+ 


1.2 




1.2 




DCO+ 


1.0 




1.0 




N2H+ * 


0.9 




1.1 




N2D+ * 


3.6 10" 


-4 


6.2 10" 


-4 


HCS+ 


1.0 




0.9 




H3O+ 


1.0 




1.5 




H2D+ * 


1.2 10" 


-4 


3.1 10" 


-4 


predictions 


X Q 


5.3 10" 


-9 


1.1 10- 


-8 


x{Et) 


1.9 10" 


-11 


5.0 10" 


-11 


x(H2D0+) 


7.8 10" 


-13 


1.2 10" 


-12 




4.5 10" 


-10 


4.2 10" 


-10 


x(NH2D) 


1.9 10" 


-13 


1.7 10" 


-13 



* Compared with the observed upper limit. 
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Table 9: Ratios of deuterated to hydrogenated ion abundances obtained in the models 



R ri t, 1 O 


NGC 2264 




W 3 IRS5 






Mod. i 


Mod. 


2 


Mod. i 


Mod. 


2 


N2D+/N2H+ 


2.6 10-3 


4.4 10" 


-3 


1.4 10-4 


2.0 10" 


-4 


H2D+/H+ 


1.7 10-3 


6.6 10" 


-4 


2.8 10-4 


2.7 10" 


-4 


DCO+/HCO+ 


3.9 10-3 


3.2 10" 


-3 


2.1 10-3 


2.1 10" 


-3 


H2DO+/H3O+ 


2.9 10-3 


1.2 10" 


-3 


1.8 10-3 


1.9 10" 


-3 


NH2D/NH3 


3.7 10-4 


1.6 10" 


-4 


4.2 10-4 


4.1 10" 


-4 



8). The calculated values of a;(N2D"'") and a;(H2D"'") are much lower than the observed 
upper limits. The electron fraction is Xe = (0.5 — 1.1) 10~^ and x(M'^) = (3 — 8) 10-^. 
M"*" is the dominant ion, providing between 60 and 70% of the electrons. The predicted 
abundances of H3 , NH3 and NH2D are one or two orders of magnitude lower than 
toward NGC 2264, due to the higher values of a;(CO) and x[0) and the lower values 
of (u and x(N2)- 

The predicted H2DO+/H3O+ ratio is close to the observed ratio DCO+/HCO+ ~ 2.5 10-3, 
but is about one order of mag nitude higher than the calculated H2D+ /R^ , N2D+ /N2H+ 
and NH2D/NH3 ratios (Table 9): these low values are mostly due to the low abun- 
dance of atomic deuterium adopted in the models required to reproduce the observed 
DCO"*" and HCO"*" abundances. 

6 Conclusions and Discussion 

We have observed submillimeter lines of several molecular ions and neutral species in 
two dense clouds, NGC 2264 and W 3 IRS5, which allow accurate determination of the 
physical conditions and abundances. Uncertainties remain due to the fact that possible 
density, temperature and abundance variations have been ignored, so that the derived 
n(H2), Tkin and abundances are averaged over the beam and along the line of sight. 
The absolute abundances are furthermore uncertain due to lingering uncertainties in 
the CO/H2 abundance ratio; the relative values should be more reliable, however. 

In both clouds the lower limit on the electron fraction derived from the sum of the 
abundances of the observed ions is Xe > 2 10-^ at densities of (1 — 2) 10^ cm-3. 
A more accurate determination of Xe has been made, using a simple steady state 
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chemical model in which the most important formation and destruction routes of the 
ions are taken into account, and in which as many species as possible are constrained by 
observations. Although there are still uncertainties in the values of the reactions rates, 
especially for dissociative recombination, we obtained agreement with the observations 
to within a factor of two. The possibility exists that a chemical steady state has not 
been reached in the observed sources, in which case comparison of observations with 
models are more difficult since the initial conditions and the age of the cloud are 
unknown. However, time-dependent models such as those of Millar et al. (f99f ) and 
Shalabiea & Greenberg (f995) show little variation in the major ion abundances and 
the electron fraction with time. 

More elaborated models using a full blown gas phase network of more than 2000 
reactions and including gas-grain interactions give comparable results: for a density 
n(H2) = 5 10^ cm~"^ and elemental gas-phase abundances close to those derived for 
NGC 2264 and W 3 IRS5, Shalabiea & Greenberg (1995) obtain x,^ ~ (1.5 - 3.5) 10"^ 
using a pure gas phase model or a dust /gas model (see their Model A). Similar results 
are presented by Le Bourlot et al. (1993) with a pure gas-phase model, as well as by 
Millar (1994, private communication). 

We did not consider in this work the interaction between electrons and grains. As 
discussed e.g. by d'Hendecourt et al. (1985) this process leads in dense clouds to grains 
with an average charge of — e, due to the inefficient photoelectric effect. Because the 
grain abundance is typically 10~^^, this results in a decrease in the gas phase electron 
fraction by only a few percent at most. Moreover, this situation may be modified if 
the internally produced ultraviolet photons due to the excitation of H2 by cosmic rays 
are taken into account, as suggested by preliminary results of Cervetto et al. (1995): 
this ultraviolet radiation could be strong enough to give rise to a photoelectric effect 
sufficient to drive the charge distribution of grains toward positive values. 

The ionization balance in dense clouds may also be modified by the presence of large 
molecules such as PAHs. In particular, the density of electrons may be reduced and 
negative PAH~ ions may contain most of the negative charge if the total PAH abun- 
dance and the thermal sticking coefficient of electrons on neutral PAH molecules are 
high (Omont 1986, Lepp & Dalgarno 1988, AUamandola et al. 1989). This means that 
the electron fractions derived in the preceding section are upper limits, but that the 
degree of ionization could be higher, typically of the order of the total fraction of PAH 
which is generally assumed to be xpah ~ 10"'' — 10~^, for PAH molecules containing 
between 1 and 10% of the total carbon and formed with about 80 C atoms. The re- 
duced electron fraction would imply an increase in the abundances of those ions which 
are predominantly destroyed by dissociative recombination - HsO"*", HCS"*", HCO"*" and 
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DCO"*" in our case - , although the (slower) destruction by mutual neutralization with 
PAH~ ions diminishes the effect. Thus, the numerical results of Section 5 could be 
slightly modified if PAH molecules are included in the modeling. However, the large 
uncertainties in the total PAH abundance in dense clouds and in the rates involving 
these species prevent clear conclusions. In particular, Mendoza-Gomez et al. (1995) 
argue that PAH molecules do not survive accretion onto grains inside dense clouds. 

In summary for dense clouds with n(H2) ~ 10^ cm~"^, the inferred electron fraction lies 
in the range 2 10~^ < Xe < 4 10~®. Ionized metals such as Mg"*" or Fe"'" could provide 
up to 60% of the electrons. These values of Xe are consistent with those commonly 
used in the literature and represented by the simple formula Xg ~ 10~^ n(H2)~^^^ 
(McKee 1989). Comparison of our results to those obtained in lower density, colder 
clouds (Schilke et al. 1991), does not show a clear dependence of Xe on density between 
10"* and 10^ cm~"^, however. In fact, the analysis of the preceeding section shows that 
other parameters which can fiuctuate from cloud to cloud also play a role. 
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